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CN ; Abstract 

G> ■ We have obtained evidence of gamma-ray emission above 1 TeV from PSR1706-44, using 

p i' a ground-based telescope of the atmospheric Cerenkov imaging type located near Woomera, 

O 
H 



South Australia. This object, a 7-ray source discovered by the COS B satellite (2CG342- 
c§ ! 02), was identified with the radio pulsar through the discovery of a 102 ms pulsed signal 
with the EGRET instrument of the Compton Gamma Ray Observatory. The flux of the 

>< ' 

\-i ' present observation above a threshold of 1 TeV is ~ 1 • 10 photons cm ^ s ^, which is 
two orders of magnitude smaller than the extrapolation from GeV energies. The analysis is 
not restricted to a search for emission modulated with the 102 ms period, and the reported 
flux is for all 7-rays from PSR1706-44, pulsed and unpulsed. The energy output in the TeV 
region corresponds to about 10^^ of the spin down energy loss rate of the neutron star. 
Subject headings: gamma rays: observations - pulsars: individual (PSR1706-44) - supernova 
remnants 
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1. Introduction 

To date, five young pulsars - the Crab (Swanenburg et al. 1981; Nolan et al. 1993), Vela 
(Swanenburg et al. 1981), Geminga (Bertsch et al. 1992), PSR1706-44 (Kniffen et al. 1992; 
Thompson et al. 1992) and PSR1055-52 (Fierro et al. 1993) - have been discovered to be 
bright in the hundreds of MeV to GeV 7-ray energy region. Inspection of these five reveals 
that the longer the pulsar period, the greater the fraction of the total spin down energy 
loss that is emitted in high energy 7-rays. The Crab nebula, containing the pulsar with the 
shortest period of these five, is so far the only established source in which 7-ray emission 
extends to the very high energy region near 1 TeV (Weekes et al. 1989; Vacanti et al. 1991) 
and even higher (Baillon et al. 1993; Tanimori et al. 1994). More observational data are 
needed to understand the processes in these objects which result in greater energy losses at 
very high energies than any other. 

The pulsar PSR1706-44 has been observed with the 3.8m diameter Cerenkov imaging 
telescope of the CANGAROO Collaboration (Patterson and Kifune 1992; Kara et al. 1993; 
Gregory et al. 1990) near Woomera, South Australia. This object was detecteded as a 7-ray 
source (2CG342-02) with the COS B satellite (Swanenburg et al. 1981) and then identi- 
fied with the radio pulsar through the discovery of 102 ms pulsed signal with the EGRET 
instrument of the Compton Gamma Ray Observatory (Kniffen ct al. 1992; Thompson et 
al. 1992). The slowing down of the pulsar period indicates that the pulsar has a charac- 
teristic age of 17,000 years. X-ray emission was also found with the ROSAT X-ray satellite 
(Becker et al. 1992), and this may suggest an X-ray synchrotron nebula associated with 
the pulsar PSR1706-44 similar to the case of the Crab. McAdam, Osborne and Parkinson 
(1993) recently suggested a possible association of the pulsar with a shell- type supernova 
remnant. Although these features are somewhat different from the Crab case, a common 
mechanism of rapid rotation of pulsar magnetosphere is likely to cause VHE 7-ray emission 
also in PSR1706-44. Cheng and Ding (1994) discussed the energy spectra of 7-ray pulsars at 
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100 MeV to 10 GeV energies to fit them with parameters based on the 'outer gap' model of 
the pulsar magnetosphere. The best fit value they have found so far for the 7-ray emission 
distance from the neutron star suggests that PSRl 706-44 should emit double pulses and 
could be a TeV 7-ray emitter. 

Evidence of TeV 7-ray emission was found in data from July and August 1992 (Ogio et 
al. 1993) and in order to confirm this result observations were undertaken in 1993 (Kifune 
et al. 1993a). This paper presents the results from these two years of observations. 

2. Instrumentation and Analysis Method 

Atmospheric Cerenkov telescopes detect optical Cerenkov photons emitted as air showers 
initiated by very high energy 7-rays and cosmic rays travel through the upper atmosphere. 
The total number of Cerenkov photons is proportional to the 7-ray energy. The 3.8 m, 
alt-azimuth mounted, telescope (Hara et al. 1993) has a ~ 1 TeV threshold for detecting 
7-rays, and is equipped with a multi-pixel camera of 220 photomultiplier tubes, each of 
which views a 0.°12 x 0.°12 area of the sky, with a total camera field of view of about 3°. 
The camera measures the number of photo-electrons and the pulse arrival time in each tube. 
The event trigger requires more than three photomultiplier tubes to have a signal of > 3 
photo-electrons and the total number of photo-electrons detected to be larger than about 
20. The arrival time measurement in each photomultiplier tube is useful for distinguishing 
the background, randomly incident, light in contrast to the nearly simultaneous Cerenkov 
photons. The night sky background hght is estimated to be about 0.05 photo electrons 
per 10 ns per photomultipher tube (Hara et al. 1993), and, thus, the background light 
which affects the Cerenkov light image consisting typically of 10 photomultiplier tubes is 
~ 1 photo-electron. A bright star, of magnitude 3.3 {r} Scorpii), is located about 1.4° from 
PSR1706-44, within the field of view of the camera. By monitoring the singles count rate in 
each phototube, we were able to monitor the passage of this star around the field of view, 
and the pointing of the telescope was calibrated. The object PSR1706-44 was kept within 
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0.°1 of the centre of the field of view. 

Monte Carlo simulations indicate that the image shape of Cerenkov hght from a 7-ray 
shower can be approximated by an ellipse, which is elongated towards the source position. 
A cosmic ray shower produces a broader, more irregularly shaped image. The image can be 
characterised by parameters such as the 'width', 'length', 'distance' (the distance between 
the centroid of the image and the centre of the field of view) and the orientation of the 
image, 'a'. The parameter a is the angle between the major axis of the elliptical image and 
the line joining the nominal source position (usually the centre of the field of view) to the 
centroid of the observed image. The definition of these parameters used in our analysis is the 
same as those of Weekes et al. (1989). The Whipple detections of very high energy 7-rays 
from the Crab nebula (Vacanti et al. 1991; Lewis et al. 1993) and the galaxy Markarian 421 
(Punch et al. 1991) have demonstrated that the 7-ray signal can be distinguished from the 
cosmic ray background using the image shape. They have shown that 7-ray events from a 
point source appear as a peak near the origin on a flat background distribution when event 
rate is plotted as a function of a. The size of this peak will be affected by the choice of the 
nominal source position used in the determination of a, so that the position of a 7-ray point 
source within the field of view can be determined to an accuracy of ~ 0.°1. 

The procedures of analysing Cerenkov images apphed to telescopes with different charac- 
teristics can be examined by using the Crab pulsar /nebula as a standard candle of very high 
energy 7-rays. The result of the Crab observations (Tanimori et al. 1994) with the present 
telescope was found to be compatible with the energy spectrum reported by the Whipple 
group (Lewis et al. 1993). 

3. Observed Data 

The image parameters were calculated for the events that are located well inside, but 
not at the centre of, the field of view, selecting the events with 'distance' smaller than 0.°9 
and greater than the 'length'. Those events with compact images (narrower than energy 
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dependent values typically of 0.°18 for 'width' and 0.°45 for 'length') were then selected to 
enrich any 7-ray component present. These values for selecting data differ slightly from those 
applied to the Crab data (0.°14 for 'width' and 0.°33 for 'length' (Tanimori et al. 1994)). 
Larger discrimination values are used in the present analysis as the image size is smaller at 
the larger zenith angles of the Crab observations (Kifune et al. 1993b). The discrimination 
values employed here result in a similar fraction of 7-rays passing the selection cut as for 
the Crab analysis (Tanimori et al. 1994). The average image size of the Crab observations 
was about 3/4 of that at the zenith, consistent with a Monte Carlo simulation for inclined 
showers (Tsukagoshi 1994). Comparison with simulations as well as the dependence of the 
chosen selection criteria on the total yield of photons will be described in a subsequent paper. 

The number of analysed events is plotted in Figure 1 as a function of the parameter 
'a'. Figure la is the plot for all the data sets combined, and Figure lb, Ic and Id are from 
the 1992 data set, July 1993 data set and August 1993 data set, respectively. In order to 
monitor the cosmic ray background in 'on-source' data, 'off-source' observation runs were 
also done. In the on-source run, PSR1706-44 was tracked at the centre of the field of view, 
whereas for off-source runs a point at the same declination as PSRl 706-44 but offset in right 
ascension was tracked. The solid line in the figure indicates the on-source data, and the 
dotted line shows the off-source data. A peak of events is seen at o; ~ 0° for the on-source 
data but not for the off-source data. In 1992, 18 hours of off-source data and 42 hours of 
on-source data were recorded and the results shown in Figure lb. For the 1993 data set, the 
observation time is 42 hours for both on- and off-source runs. The event rates of on- and 
off-source observations were equal to each other at a > 30°. The peak at a ~ 0° is thus 
due to an excess of gamma ray events above the background of isotropic cosmic rays. The 
excess at a ~ 0° is present also when no selection cut was made on 'width', and becomes 
more prominent as narrower widths arc selected, as expected if a 7-ray signal is present. The 
width of the observed peak is consistent with that estimated from Monte Carlo simulations 
for 7-ray emission from a point source. 
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An analysis was performed in which the 'source' position was dehberately shifted (in 
software) from PSR1706-44. This procedure leads to a change in calculated values of 'a'. 
We then studied how the peak strength at a ~ 0° for these 'false sources' varied as a function 
of the shift. The contour map of the significance of the peak (calculated from the number 
of events in a < 10° compared with a > 30°) is shown in Figure 2 as a function of the 
position around PSR1706-44. Since the telescope is alt-azimuth mounted, the field of view 
rotates around its center with time. The 'de-rotation' to correct for this effect was made to 
the data in software. The highest significance is obtained when the source location is set at 
the centre of the field of view, the true position of PSR1706-44. 

The number of data samples in each map in Figure 2 is 1681. We would therefore expect a 
few ±3(7 statistical fluctuations. However, deviations of about Aa are evident near the outer 
regions of the map. In the analysis to obtain Figure 2, the requirement that 'distance' be 
smaller than 0.°9 was applied to the 'distance' of the centroid of image from each 'false source' 
position as well as the usual 'distance' from the center of field of view. This procedure resulted 
in unequal areas in the field of view which contribute to the data for different assumed source 
positions, and led to poorer statistics for the 'false source' positions near the edge of the field 
of view. A strict comparison between 'source' and 'false source' should be made to similarly 
sized data-bases (as proposed in Fomin et al. 1994). In the present case, this condition is 
met only at the center of the field of view. However, the inequality of the angular area is 
not large for the positions near the center, and 'false sources' in this central region do not 
produce fluctuations larger than expected from statistics. The significance obtained for the 
true position of PSRl 706-44 is the only one that exceeds statistical expectations near the 
center. The existence of 4(7 fiuctuations near the edges of the map suggests a limitation to 
this contour map method when applied to 'false source' locations considerably removed from 
the center of the field of view. 

The selection of 'distance' less than 0.°9 is found effective in suppressing any effects of the 
tendency for images located near the edge of field of view to appear elongated along the outer 
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boundary of the camera, i.e. having values of a near 90°. Enrichment of a > 30° events as 
a consequence of this effect would result in the number of events at ct ~ 0° becoming smaller 
than the rate estimated from a of larger values. 

We examined the data for the presence of any systematic effect. The events that con- 
tribute to the peak are distributed over the full range of 7-ray energy as estimated from 
the total yield of Cerenkov photons. Thus, the peak at a; ~ 0° is not caused by a spurious 
effect near the hmit of detection sensitivity. The plot of the centroid positions in the image 
plane shows a uniform distribution indicating no correlation with any specific position, such 
as that of rjScorpii. This star is offset from PSR1706-44 by 0.°64 in right ascension and by 
l.°25 in declination: a location beyond the upper right edge in the map of Figure 2. The 
separation of this star from PSR1706-44 is greater than the 0.°9 discrimination value on 
'distance', and, thus, it is unlikely that the analysed images are affected by rjScorpii. 

The event rate of the cosmic ray background has a flat distribution as a function of the 
parameter a. By estimating the background from the rate in the range a — 30° ~ 90°, the 
excess counts within the bins with a < 10° in Fig. la has a significance of 12 a. The number 
of events in the peak is constant within statistical fluctuations when the data arc divided into 
shorter sub-sets. The 7-ray intensity is, therefore, consistent with steady emission during 
the observation period. 

The effective detection area is a product of the area of the Cerenkov light pool and the 
detection efficiency of the telescope. The detection efficiency is a function of 7-ray energy and 
(the location of telescope in) the lateral spread of a Cerenkov hght pool, and also depends on 
the trigger conditions and on the selection criteria of events for analysis. Results of a Monte 
Carlo simulation study under way to estimate the detection efficiency for these observations, 
suggest an effective area of approximately 5x 10^ cm^ at a 7-ray energy of 1 TeV. Adopting 
this area, the flux of 7-ray events is 7 • rj'^ x 10~^^ cm~^ s~^, where r] is the efficiency for 
7-rays passing through the data reduction procedures. The number of excess counts which 
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constitute the peak at a ~ 0° was compared between the two cases with and without the 
selection cuts on 'width' and 'length'. Ninety per cent of the peak events survive after the 
cuts, similar to the case of the Crab data (Tanimori et al. 1994). By putting r] equal to 0.9, 
we obtain a 7-ray flux of 8 x 10~^^ cm~^ s~^, integrated above the threshold energy of 1 TeV. 
We estimate that the systematic error in the flux and 7-ray energy may be as large as ~ 30 
% from uncertainties in the adopted values of effective area and efficiency rj. Searches for 
pulsed emission correlated with the 102 ms period (in preparation for publication) indicate 
that it is unlikely that a major portion of the detected flux is pulsed. 

4. Discussions 

The present yield of 7-rays corresponds to ~ 1 • 10~^^ erg cm~^ s~^ at 1 TeV. For 
a distance of 1.5 kpc to PSR1706-44, the TeV luminosity is ~ 3 ■ 10^^ erg s~^ assuming 
isotropic emission. This is 10~^ of the loss rate of the total spin down energy of 3.4 • 10^^ 
erg s~^ of the neutron star. The TeV luminosity is smaller by an order of magnitude than 
the luminosity of 2.6 ■ 10^^ erg s~^ in the 100 MeV - 10 GeV range, but greater than the 
X-ray luminosity of 1 • 10^^ erg s"^ in the 0.1 - 2.4 keV range detected with the ROSAT. 
The energy spectrum of the pulsed signal reported in Thompson et al. 1992 is very flat up 
to several GeV with a differential spectral index of —1.7. When extrapolated to 1 TeV the 
integral flux is as high as 1.6 x 10~^ photons cm"^ s~^ (time-averaged over the period). The 
flux at 1 TeV from the present observation is smaller by two orders of magnitude than this 
extrapolated value. 

We note that Nel et al. (1993) have reported an upper hmit of 5.8 x 10~^^ cm~^ s~^ 
for > 2.6 TeV on the 7-ray flux that is modulated with the pulsar period. However, our 
analysis is not restricted to a search for pulsed emission, and is sensitive to any 7-rays from 
PSR1706-44 - pulsed and unpulsed. 

It has been reported by McAdam, Osborne and Parkinson (1993) that PSR1706-44 is 
located near the edge of a shell-shaped supernova remnant of about a half degree diameter. If 
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the production mechanism of TeV 7-rays is similar to that of the Crab, then pulsar PSR1706- 
44 may also have a synchrotron X-ray nebula like the Crab, and the ROSAT result may be 
an indication of this. Although the peak of the a-distribution is consistent with a point 
source, the angular resolution of the telescope, ~ 0.1°, corresponding to a size <, 3 pc at 
a distance of 1.5 kpc, does not preclude such extended emission. Multi- wavelength studies 
of this object will be needed in order to investigate the production mechanism for the very 
high energy 7-rays reported here. 
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Figure legends 

FIG. 1 The number of observed events as a function of the orientation angle a of the 
Cerenkov light image. On- and off-source data are indicated by the solid and dotted hues, 
respectively. (The on-source and off-source observing times are not equal.) (a) All data, (b) 
The data from 1992. (c) The data from July 1993. (d) The data from August 1993. 

FIG. 2 The contour map of statistical significance for various directions in the sky 
around the PSR1706-44 position. In calculating the significance, the assumed position of 
the '7-ray point source' was artificially varied around the known position of PSR1706-44 and 
the significance of any resulting peak in the a distribution calculated. The result shown in 
the figure is from the data set of July and August 1993. Contours of an equal significance 
are shown in the lower figures against the two dimensional directions in units of degrees. 
The upper figures are the corresponding lego plots. In order to have the same scale both 
in on- and off-source plots to facilitate comparison, an artificial highest significance of 8.5 a 
was inserted in the farthest left-hand bin. The figures on the left are for on-source data and 
those on the right for off-source data. 
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